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Abstract. A new era in fundamental physics began when pulsars were discovered in 
1967. Soon it became clear that pulsars were useful tools for a wide variety of physical 
and astrophysical problems. Further applications became possible with the discovery 
of the first binary pulsar in 1974 and the discovery of millisecond pulsars in 1982. Ever 
since pulsars have been used as precise cosmic clocks, taking us beyond the weak-field 
limit of the solar-system in the study of theories of gravity. Their contribution is crucial 
as no test can be considered to be complete without probing the strong-field realm of 
gravitational physics by finding and timing pulsars. This is particularly highlighted by 
the discovery of the first double pulsar system in 2003. In this review, I will explain 
some of the most important applications of millisecond pulsar clocks in the study of 
gravity and fundamental constants. 



1 Introduction 

The title of this volume, "Astrophysics, Clocks, and Fundamental Constants", 
would also be a suitable title for this contribution describing the use of radio 
pulsars in the study of fundamental physics. Indeed, pulsar astronomy is an 
extraordinary discipline which removes the distinction between physics and as- 
trophysics that is often made. Such a distinction may be justified by the fact 
that in a terrestrial laboratory we can modify the experimental set-up and con- 
trol the environment. In contrast, in astrophysical experiments we remain an 
observer, deriving all our information simply from observing photons and their 
properties. Thereby, terrestrial experiments are typically more precise and, most 
importantly, can be reproduced - at least in principle - in any other laboratory 
on Earth. However, when probing the limits of our understanding of fundamental 
physics, we often have to study conditions that are too extreme to be encoun- 
tered on Earth. One may take the experiment into space, like "LISA" , "STEP" 
or "Gravity Probe-B" , but even then we are limited, particularly if we want to 
study gravity. While solar system tests provide a number of very stringent tests 
for general relativity, none of the experiments made or proposed for the future 
will ever be able to test the strong field limit. For such studies, pulsars are and 
will remain the only way to test and enhance our understanding. Additionally, 
pulsars not only provide us with the only means to perform strong-field experi- 
ments, but these experiments are also amazingly precise. It is this unique aspect 
that I will review in the following. The interested reader may also consult the 
excellent reviews by Will and Turyshev |2] on PPN formalism, by Wex |3| 
and Stairs 0] on strong gravity tests, and by Lorimer j^j on pulsars in general. 
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2 Pulsars 

Pulsars are highly magnetized, rotating neutron stars which emit a narrow radio 
beam along the magnetic dipolc axis. As the magnetic axis is inclined to the 
rotation axis, the pulsar acts like a cosmic light-house emitting a radio pulse 
that can be detected once per rotation period when the beam is directed towards 
Earth (Fig. 0). For some very fast rotating pulsars, the so-called millisecond 
pulsars, the stability of the pulse period is similar to that achieved by the best 
terrestrial atomic clocks. This is not surprising if we consider that they have 
large rotational energies of E = 10 43-45 J and low energy loss rates. Using these 
astrophysical clocks by accurately measuring the arrival times of their pulses, a 
wide range of experiments is possible, some of which are presented here. While 
it is not of utmost importance for the remainder of this review how the radio 
pulses are actually created, we will consider some of the basic pulsar properties 
below. 

2.1 Pulsars as Neutron Stars 

Pulsars are born in supernova explosions of massive stars. Created in the col- 
lapse of the stars' core, neutron stars are the most compact objects next to black 
holes. From timing measurements of binary pulsars (see Section I5.2JI . we deter- 
mine the masses of pulsars to be within a narrow range of (1.35 ± 0.04) M |S]. 
Modern calculations for different equations of state produce results for the size 
of a neutron star quite similar to the very first calculations by Oppenheimer & 
Volkov 0, i.e. about 20 km in diameter. Such sizes are consistent with indepen- 
dent estimates derived from modelling light-curves and luminosities of pulsars 
observed in X-rays (e.g. [E])- 

As rotating magnets, pulsars emit magnetic dipole radiation as the dominant 
effect for an increase in rotation period, P, described by P. Equating the cor- 
responding energy output of the dipole to the loss rate in rotational energy, we 
obtain an estimate for the magnetic field strength at the pulsar surface from 

B s = 3.2 x 1Q 19 VP? Gauss, (1) 

with P measured in s and P in s s _1 . Sometimes twice the value is quoted to 
reflect the field at the poles. Typical values arc of order 10 12 G, although field 
strengths up to 10 14 have been observed 0. Millisecond pulsars have lower field 
strengths of the order of 10 8 to 10 10 Gauss which appear to be a result of their 
evolutionary history (see Section [3J. These magnetic fields are consistent with 
values derived from X-ray spectra of neutron stars where we observe cyclotron 
lines HOI. 

2.2 Pulsars as Radio Sources 

The radio signal of a pulsar is usually weak, both because the pulsar is distant 
and the size of the actual emission region is small. Estimates range down to a few 
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Fig. 1. (left) A pulsar is a rotating, highly magnetised neutron star. A radio beam 
centred on the magnetic axis is created at some height above the surface. The tilt 
between the rotation and magnetic axes makes the pulsar in effect a cosmic lighthouse 
when the beam sweeps around in space, (right) Individual pulses vary in shapes and 
strength (top) average profiles are stable (bottom). The typical pulse width is only 
~4% of the period. 



metres, resulting in brightness temperatures of up to 10 37 K Such values 
require a coherent emission mechanism which, despite 35 years of intensive re- 
search, is still unidentified. However, we seem to have some basic understanding, 
in which the magnetized rotating neutron star induces an electric quadrupolc 
field which is strong enough to pull out charges from the stellar surface (the 
electrical force exceeds the gravitational force by a factor of ~ 10 12 !). The mag- 
netic field forces the resulting dense plasma to co-rotate with the pulsar. This 
magnetosphere can only extend up to a distance where the co-rotation veloc- 
ity reaches the speed of light 1 . This distance defines the so-called light cylinder 
which separates the magnetic field lines into two distinct groups, i.e. open and 
closed field lines. The plasma on the closed field lines is trapped and co-rotates 
with the pulsar forever. In contrast, plasma on the open field lines can reach 
highly relativistic velocities and can leave the magnetosphere, creating the ob- 
served radio beam at a distance of a few tens to hundreds of km above the pulsar 
surface (e.g. see Fig. 0). 



1 Strictly speaking, the Alfven velocity will determine the co-rotational properties of 
the magnetosphere. 
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Most pulsars are not strong enough for us to allow studies of their individual 
radio pulses. Then, only an integrated pulse shape, the "pulse profile" , can be 
observed. Individual pulses reflect the instantaneous plasma processes in the 
pulsar magnetosphere, resulting in often seemingly random pulses (see Fig. EJ. 
In contrast, the average pulse profile reflects the global constraints mostly given 
by a conal beam structure and geometrical factors and is thereby stable. It is 
this profile stability which allows us to time pulsar to high precision. 
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Fig. 2. The P — P-diagram for the known pulsar population. Lines of constant char- 
acteristic age and surface magnetic field are shown. Binary pulsars are marked by a 
circle. The solid line represents the pulsar "death line" enclosing the "pulsar graveyard" 
where pulsars are expected to switch off radio emission. 



3 A Pulsar's Life 

The evolution in pulsar period, P, and slow-down, P, can be used to describe the 
life of a pulsar. This is usually done in a (logarithmic) P-P-diagram as shown 
in Fig. [21 where we can draw lines of constant magnetic field (see Eq. JIJ) an d 
constant "characteristic age" estimated from 

P v 

t=^ = , 2) 

2P 2v V ' 

using either period, P, or the spin frequency, v, and their derivatives in standard 
units. This quantity is a valid estimate for the true age under the assumption 
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that the initial spin period is much smaller than the present period and that the 
spin-down is fully determined by magnetic dipole braking. While it had been 
assumed that pulsars are born with birth periods similar to that estimated for 
the Crab pulsar, Pq = 19 ms [13] . recent estimates for a growing number of 
pulsars suggest a wide range of initial spin periods from 14 ms up to 140 ms 
[Tl] . Pulsars are therefore born in the upper left area of Fig. |3 and move into 
the central part where they spent most of their lifetime. 

3.1 Normal Pulsars 

Most known pulsars have spin periods between 0.1 and 1.0 sec with period deriva- 
tives of typically P = 10~ 15 s s _1 . Selection effects are only partly responsible for 
the limited number of pulsars known with very long periods, the longest known 
period being 8.5 s |15|. The dominant effect is due to the "death" of pulsars 
when their slow-down has reached a critical state. This state seems to depend 
on a combination of P and P which can be represented in the P — P-diagram 
as a "pulsar death-line" . To the right and below this line (see Fig. the elec- 
tric potential above the polar cap may not be sufficient to produce the particle 
plasma that is responsible for the observed radio emission. While this model can 
indeed explain the lack of pulsars beyond the death-line, the truth may be more 
complicated as the position of the 8.5-sec pulsar deep in the "pulsar graveyard" 
indicates. Nevertheless, it is clear that the normal life of radio pulsars is limited 
and that they die eventually after tens to a hundred million years. 

3.2 Millisecond Pulsars 

Inspecting the approx. 1600 sources shown in the P — P-diagram, it is obvious 
that the position of a sub-set of about 100 pulsars located in the lower left 
part of the diagram cannot be explained by the above picture of normal pulsar 
life. Instead, these pulsars simultaneously have small periods (of the order of 
milliseconds) and small period derivatives, P <; 10~ 18 s s _1 . They appear much 
older than ordinary pulsars (see Eq. (0) and, indeed, these so-called "millisecond 
pulsars" represent the oldest population of pulsars with ages up to ~ 10 10 yr. 
A model for their evolutionary history was proposed soon after the discovery of 
PSR B1937+21 by Backer et al. in 1982 [B . This first millisecond pulsar has a 
period of only 1.56 ms and remains the pulsar with the shortest period known. 

It is suggested that millisecond periods are obtained when mass and thereby 
angular momentum is transferred from an evolving binary companion while it 
overflows its Roche lobe In this model, millisecond pulsars are recycled 
from a dead binary pulsar via an X-ray binary phase. This model implies a 
number of observational consequences: a) most normal pulsars do not develop 
into a millisecond pulsar as they have long lost a possible companion during 
their violent birth event; b) for surviving binary systems, X-ray binary pulsars 
represent the progenitor systems for millisecond pulsars; c) the final spin period 
of recycled pulsars depends on the mass of the binary companion. A more massive 
companion evolves faster, limiting the duration of the accretion process; d) the 
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majority of millisecond pulsars have low-mass white-dwarf companions as the 
remnant of the binary star. These systems evolve from low-mass X-ray binary 
systems (LMXBs); e) high-mass X-ray binary systems (HMXBs) represent the 
progenitors for double neutron star systems (DNSs). DNSs are rare since these 
systems need to survive a second supernova explosion. The resulting millisecond 
pulsar is only mildly recycled with a period of tens of millisecond. 

The properties of millisecond pulsars and X-ray binaries are consistent with 
the described picture. For instance, it is striking that ~ 80% of all millisecond 
pulsars are in a binary orbit while this is true for only less than 1% of the 
non-recycled population. For millisecond pulsars with a low-mass white dwarf 
companion the orbit is nearly circular. In case of double neutron star systems, the 
orbit is affected by the unpredictable nature of the kick imparted onto the newly 
born neutron star in the asymmetric supernova explosion of the companion. If 
the system survives, the result is typically an eccentric orbit with an orbital 
period of a few hours. As we will see, both types of system can be used to test 
different aspects of gravitational theories. 

4 Pulsars as Clocks 

By measuring the arrival time of the received ticks of the pulsar clock very 
precisely, we can study effects that determine the propagation of the pulses in 
four-dimensional space-time. Millisecond pulsars are the most useful objects for 
these investigations: their pulse arrival times can be measured much more accu- 
rately than for normal pulsars (the measurement precision scales essentially with 
spin frequency) and their rotation is much smoother, making them intrinsically 
better clocks. Specifically, they do not exhibit rotational instabilities known for 
normal pulsars, namely "timing noise" and "glitches" . Glitches are associated 
with young pulsars and they represent a sudden increase in rotation frequency 
that is probably caused by an abrupt change in the internal structure of the 
neutron star. The origin of timing noise is much less understood. It manifests 
itself in a quasi-random walk in one or more of the rotational parameters on 
timescales of months to years. Again, it appears mostly for young pulsars and 
scales with some power of the period derivative, P. Hence, millisecond pulsars 
generally do not show timing noise, although it has been detected for few sources 
such as PSR B1937+21 albeit on a much smaller amplitude scale than for 
normal pulsars. 

4.1 Time Transfer 

In order to study effects that change the pulse travel time, we first have to find 
an expression that describes the pulsar rotation in a reference frame co-moving 
with the pulsar. We start by expressing the spin frequency of the pulsar in a 
Taylor expansion, 



v(t) = v + i>o(t- *o) + ^o{t - t ) 2 + 



(3) 
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where vq is the spin frequency at reference time to, i.e. Vq — v(to) = 1/Pq with 
Po being the corresponding pulse period. While vq and its derivatives refer to 
values measured at a certain epoch, i>, v are determined by the physical process 
responsible for the pulsar slow-down and should, in principle, be constant for 
most time-spans considered. We expect a relation 

v = -const. v n ', (4) 

and hence i) — —const, x n x v n i> where the "braking index", n, has a 
value of n = 3 for magnetic dipole braking, relating to Eq. Measuring v 
can yield the braking index via n — vv/v 2 so that the assumption of dipole 
braking can be tested. Timing noise can mimic a significant but time-varying 
value of v that reflects timing noise rather than regular spin-down. In these 
cases, derived braking indices are meaningless in terms of global spin-down. For 
most millisecond pulsars v is too small to be of significance although some source 
show a non-zero v due to timing noise. 

Relating the spin frequency to the pulse number N, we find 

N = N + v {t - t ) + |z> (t ~ t ) 2 + -i> {t - to) 3 + ■■■ (5) 

where Nq is the pulse number at the reference epoch to. If to coincides with the 
arrival of a pulse and the pulsar spin-down is accurately known, the pulses should 
therefore appear at integer values of N when observed in an inertial reference 
frame. 

Our observing frame is not inertial, as we are using telescopes that are located 
on a rotating Earth orbiting the Sun. Before analysing corresponding TOAs 
measured with the observatory clock ( "topocentric arrival times"), we need to 
transfer them to the centre of mass of the solar system as the best approximation 
to an inertial frame available. By using such "barycentric arrival times" , we can 
easily combine transferred topocentric TOAs measured at different observatories 
at different times. The transformation of a topocentric TOA to a barycentric 
arrival time, tssB, is given by 

tsSB = ^topo — ^0 + icorr — D/ f 2 , 

+ ^Rocmor,0 + ^Shapiro,© + ^Einstein,© , 
~t~ -^RoemerjBin "1" ^Shapiro, Bin "T" ^Einstein, Bin 

We have split the transformation into three lines. The first two lines apply to 
every pulsar whilst the third line is only applicable to binary pulsars. We discuss 
each term in detail. 

Clock and frequency corrections The observatory time is typically main- 
tained by local H-maser clocks that are compared to UTC (NIST) by the Global 
Positioning System (GPS). Offsets are monitored and retroactively applied as 
clock corrections, t COII , in the off-line analysis which often uses UTC of BIPM 
as time standard. Further corrections take into account that the Earth is not 



(6a) 
(6b) 
(6c) 
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rotating uniformly, so that leap seconds are occasionally inserted into UTC to 
keep it close mean solar time. All leap second are removed from the used UTC 
time standard to produce a TOA measured in International Atomic Time (TAI) . 
The TAI is maintained as an average of a large number of selected atomic clocks 
by the Bureau International des Poids et Mesures (BIPM) , which also publishes 
a retroactive uniform atomic time standard known as Terrestrial Time, TT (for- 
merly known as Terrestrial Dynamical Time, TDT). The unit of TT is the SI 
second and may be regarded as the time that would be kept by an ideal atomic 
clock on the geoid with TT=TAI+32.184 seconds, where the offset of about 32 s 
stems from historic reasons. This time scale should be used in the final analysis 
by correcting the initially measured TO As to TT(BIPM). 

As the pulses are delayed due to dispersion in the interstellar medium, the ar- 
rival time depends on the observing frequency, /. The TOA is therefore corrected 
for a pulse arrival at an infinitely high frequency, thereby removing dispersion 
from the data (last term in Ea. lfiajl . The corresponding Dispersion Measure (DM) 
is determined during the discovery of the pulsar and can be measured accurately 
by observations at multiple frequencies. For some pulsars, the dispersion measure 
is observed to change with time. In order to avoid time- varying drifts introduced 
into the TOAs in such cases, the above term needs to be modified to include 
time-derivatives of DM, i.e. DM, DM and so on. These can be determined if 
monitoring observations at two or more frequencies are available when they pro- 
vide an estimate for the change in electron density along the line-of-sight as a 
result of 'interstellar weather'. For high-precision timing of millisecond pulsars 
such multi-frequency observations are essential. 

Barycentric corrections The terms in Eq. Ijfibjl describe the corrections nec- 
essary to transfer topocentric to barycentric TOAs. 

The Roemer delay, ZiR oemerj 0, is the classical light-travel time between the 
phase centre of the telescope and the solar system barycentre (SSB). Given a 
unit vector, s, pointing from the SSB to the position of the pulsar and the vector 
connecting the SSB to the observatory, r, we find: 

^Rocmor.0 = -- r • S - (r SSB + r EO ) ■ §. (7) 

c c 

Here c is the speed of light and we have split r into two parts. The vector, 
r ssB, points from the SSB to centre of the Earth (geocentre). Computation of 
this vector requires accurate knowledge of the locations of all major bodies in 
the Solar system and uses solar system ephemerides such as the 'DE200' or 
'DE405' published by Jet Propulsion Laboratory (JPL) JHj- The second vector 
rEO j connects the geocentre with the phase centre of the telescope. In order to 
compute this vector accurately, the non-uniform rotation of the Earth has to 
be taken into account, so that the correct relative position of the observatory 
is derived. This is achieved using appropriate UT1 corrections published by 
International Earth Rotation Service (IERS). 

The Shapiro delay, ^Shapiro,© > is a relativistic correction that corrects for 
extra delays due to the curvature of space-time caused by the presence of masses 
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in the solar system The delays are largest for a signal passing the Sun's limb 
(~ 120 /is) while Jupiter can contribute as much as 200 ns. In principle one has 
to sum over all bodies in the solar system, yielding 



where Mi is the mass of body i, rf is the pulsar position relative to it, and rf is 
the telescope position relative to that body at the time of closest approach of the 
photon (see [2]). The parameter 7 is one of the Parameterised-Post- Newtonian 
(PPN) parameters that will be discussed in Section It describes how much 
space-curvature is produced by unit rest mass and takes the value 7 = 1 in 
general relativity 2 In practice, 7 is adopted as unity and only the Sun, and in 
some cases Jupiter, are accounted for in this calculation. 

The last term in Eq. i|6bj) . /^Einstein,©; is called Einstein delay and it de- 
scribes the combined effect of time dilation due to the motion of the Earth and 
gravitational redshift caused by the other bodies in the Solar system. This time 
varying effect takes into account the variation of an atomic clock on Earth in the 
changing gravitational potential as the Earth follows its elliptical orbit around 
the Sun. The delay amounts to an integral of the expression \H\ 



where the sum is again over all bodies in the solar system but this time excluding 
Earth. The distance rf is again the distance between Earth and body i, while 
ve is the velocity of the Earth relative to the Sun. 

Relative Motion Equation I(jal6cl is sufficient to measure the clock rate as 
produced by the pulsar if no further motion or acceleration between pulsar and 
SSB occurs. If the pulsar is moving relative to the SSB, only the transverse 
component of the velocity, v t , can be observed from timing. A radial motion 
is not measurable practically (though theoretically possible), leaving resulting 
Doppler corrections to observed periods, masses etc. undetermined. The situation 
changes if the pulsar has an optically detectable companion such as a white dwarf 
for which Doppler shifts can be measured from optical spectra. In contrast, a 
transverse motion will change the vector s in Eq. Q, adding a linear time- 
dependent term to our transfer equation, and can therefore be measured as 
proper motion, \i. 

Another effect arising from a transverse motion is the Shklovskii effect, also 
known in classical astronomy as "secular acceleration". With the pulsar mo- 
tion, the projected distance of the pulsar to the SSB is increasing, leading to a 

2 Usually, this parameter is denoted without the "hat" simply as 7 (see but it is 
common practice in the study of binary pulsars to use the symbol 7 to describe the 
amount of time dilation and gravitational redshift caused by a pulsar companion. 
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correction that is quadratic in time |21| , 

do) 

Since this delay scales with inverse of the distance d to the pulsar, the correction 
is usually too small to be considered. However, it has the effect that any observed 
change in a periodicity (i.e. change in pulse or orbital period) is increased over 
the intrinsic value by 

P c d K ' 

For millisecond pulsars where P is small, a significant fraction of the observed 
change in period can be due to the Shklovskii effect. This effect also needs to be 
considered, when studying the decay of an orbital period due to gravitational 
wave emission, where the observed value is increased by the Shklovskii term. 

Similarly, any line-of-sight acceleration a of the pulsar due to an external 
gravitational field changes the observed period derivative by P/P = a/c. This 
effect is commonly observed for pulsars in globular clusters where the accelera- 
tion through the cluster's gravitational field towards our line of sight can often 
be large enough to reverse the sign of P. As a result, the pulsars appear to 
be spinning up rather than down! Such pulsars place useful constraints on the 
cluster mass distributions and the intracluster medium |22| . 

Finally, a related term that needs to be considered for nearby pulsars de- 
scribes an annual parallax given by |21j 

^ = -^(-.) 2 = ^((r.^-H 2 )- d2) 

In comparison to the more familiar positional parallax, this timing parallax cor- 
responds to measuring the time delay due to the curvature of the emitted wave- 
fronts at different positions of the Earth orbit. This effect imposes a signal with 
an amplitude of l 2 cos/3/(2cd) where I is the Earth-Sun distance and (3 is the 
ecliptic latitude of the pulsar. For a pulsar at d = 1 kpc, this delay amounts 
to only <; 1.2//S, and hence it is only measurable for a few millisecond pulsars 
where it provides a precise distance estimate. Similarly difficult to measure is 
the annual- orbital parallax for binary pulsars which manifests itself as a periodic 
change in the observed projected semi-major axis of the pulsar's orbit due to 
viewing the system from slightly different directions during the Earth's orbit. In 
contrast, a secular change of the semi-major axis due to a proper motion of the 
system on the sky has been measured for a number of binary millisecond pulsars 
(see E|). 



4.2 Pulsar Timing 

The transfer equation Eqs. (|6a|l -ll^c | contains a number of parameters which 
are not known a priori (or only with limited precision after the discovery of a 
pulsar) and need to be determined precisely in a least-squares fit analysis of the 
measured TOAs. These parameters can be categorised into three groups: 
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Fig. 3. Timing residuals for the 1.19-s pulsar B1133+16. A fit of a perfect timing model 
should result in randomly distributed residuals, shown in a). A parabolic increase in 
the residuals in part b) is obtained if P is underestimated, here by 4%. An offset in 
position produces sinusoid residuals shown in part c) where the declination has an 
error of 1 arcmin. Part d) demonstrates the effect of neglected proper motion, here of 
fi = 380 mas/yr. Note the different scales on the y-axes. 



Astrometric parameters: The astrometric parameters include the position of 
the pulsar, and its proper motion and parallax. While the position is only 
known within a telescope beam after the discovery, the precision can be 
greatly improved by timing the pulsar for about a year (full Earth orbit). 
Proper motion and parallax only become evident after a longer time-span. 

Spin parameters: These include the rotation frequency of the pulsar, v 1 and 
its derivatives J3J). 

Binary parameters: For pulsars in a binary orbit, the initial observations will 
typically show a periodic variation in observed pulse period. Five Keplerian 
parameters then need to be determined: orbital period, Pf,; the projected 
semi- major axis of the orbit, x = asini where i is the (usually unknown) 
inclination angle of the orbit; the orbital eccentricity, e; the longitude of 
periastron, ui\ and the time of periastron passage, Trj. For a number of bi- 
nary systems this Newtonian description of the orbit is not sufficient and 
relativistic corrections need to be applied, e.g. uo is replaced by u> + tot. The 
measurement of the Post- Keplerian (PK) parameters such as Co allows a com- 
parison with values expected in the framework of specific theories of gravity. 
We discuss these aspects further below. 

Given a minimal set of starting parameters, a least-squares fit is needed to match 
the measured arrival times to pulse numbers according to Eq. (JSJ. The aim is 
to obtain a phase-coherent solution that accounts for every single rotation of 
the pulsar between two observations. One starts off with a small set of TO As 
that were obtained so closely in time, that the accumulated uncertainties in the 
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starting parameters do not exceed one pulse period. Gradually, the data set is 
expanded, maintaining coherence in phase. When successful, post-fit residuals 
expressed in pulse phase show a random distribution around zero (see Fig EJ). 
After starting with fits for only period and pulse reference phase over some hours 
and days, longer time spans slowly require fits for parameters like spin- frequency 
derivative (s) and position. Incorrect or incomplete timing models cause system- 
atic structures in the post-fit residuals identifying the parameter that needs to be 
included or adjusted (see Fig.^J). The precision of the parameters improves with 
length of the data span and the frequency of observation, but also with orbital 
coverage in the case of binary pulsars. Sufficient data sets then enable measure- 
ments with amazing precision, e.g. the period determined for PSR B1937+21 is 
known to a relative precision of 10 -15 . 

Obviously, the above description shows that the process of pulsar timing 
is elaborate. Fortunately, sophisticated software packages have been developed 
that combine time transfer and the least-squares fit of the timing model. The 
three major packages are TEMPO (ATNF/Princcton University) [H], PSRTIME 
(Jodrell Bank Observatory) |2U, and TIMAPR (Pushchino Observatory/MPIfR) 
The most widely used package is TEMPO. 

5 Applications of Pulsars 

Pulsars are unique and versatile objects which can be used to study an extremely 
wide range of physical and astrophysical problems. Beside testing theories of 
gravity one can study the Galaxy and the interstellar medium, stars, binary 
systems and their evolution, solid state physics and the interior of neutron stars. 
Investigating the radio emission of pulsars provides insight into plasma physics 
under extreme conditions. In the following we will concentrate on the application 
of pulsars as clocks, paying in particular attention to tests of theories of gravity. 
Some of these tests involve studies of PPN parameters and possible related time 
variation in the Gravitational Constant, G. 

5.1 PPN parameters 

Metric theories of gravity assume (i) the existence of a symmetric metric, (ii) 
that all test bodies follow geodesies of the metric and (Hi) that in local Lorentz 
frames the non-gravitational laws of physics are those of special relativity. Under 
these conditions we can study metric theories with the Parameterised Post- 
Newtonian (PPN) formalism by describing deviations from simple Newtonian 
physics in the slow-motion and weak-field limit. This is possible in a theory- 
independent fashion, such that the only differences in these theories occur in 
the numerical coefficients that appear in the metric, characterised by a set of 
10 real- valued PPN-Parameters |2fi) . Each of the parameters can be associated 
with a specific physical effect, like the violation of conservation of momentum 
or equivalence principles, and certain values are assigned to them in a given 
theory. Thereby, comparing measured PPN parameters to their theoretical values 
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can single out wrong theories in a purely experimental way. A more complete 
description of the PPN formalism and the physical meaning of PPN parameters 
is presented by Turyshev in this volume • A recent review was given by Will pQ . 
Here we summarize the studies of those PPN parameters that can be constrained 
by pulsars. A more detailed account of related pulsar tests is given by Wex \'2713\ 
and Stairs 0]. 

Violations of the Strong-Equivalence-Principle (SEP) The Strong Equiv- 
alence Principle (SEP) is completely embodied into general relativity, while al- 
ternative theories of gravity predict a violation of some or all aspects of SEP. The 
SEP is, according to its name, stronger than both the Weak Equivalence Princi- 
ple (WEP) and the Einstein Equivalence Principle (EEP). The WEP states that 
all test bodies in an external gravitational field experience the same acceleration 
regardless of the mass and composition. While the WEP is included in all metric 
theories of gravity, the EEP goes one step further and also postulates Lorentz- 
invariance and positional invariance. Lorentz-invariance means that no preferred 
frame exists, so the outcome of a local non-gravitational experiment is indepen- 
dent from the velocity of the apparatus, while positional invariance renders it 
unimportant where this experiment is being performed. The SEP includes both 
the WEP and the EEP, but postulates them also for gravitational experiments. 
As a consequence, both Lorentz- and positional invariance should be indepen- 
dent of the gravitational self-energy of the bodies in the experiment. Obviously, 
all bodies involved in terrestrial lab experiments possess only a negligible frac- 
tion of gravitational self-energy, so that tests of SEP require the involvement of 
astronomical objects. 

A violation of SEP means that there is a difference between gravitational 
mass, M g , and inertial mass, Mi. The difference can be written as 

MjL = l + 5(e) = l + V e + 0(e 2 ), (13) 

where e is the gravitation self-energy in units of rac 2 and r\ is a parameter 
characterising the violation of SEP. The latter parameter was introduced by 
Nordvedt (1968) who suggested studying the Earth-Moon system to test for vi- 
olations of SEP. Due to their different self-energy (Earth: e ~ —4.6 x 10~ 10 , 
Moon: e ~ —0.2 x 1CT 10 ), Earth and Moon would fall differently in the exter- 
nal gravitational field of the Sun, leading to a polarization of the Earth-Moon 
orbit ( "Nordvedt-effect" ) . Lunar-laser-ranging experiments can be used to put 
tight limits on r\ which is a linear combination of PPN parameters represent- 
ing effects due to preferred locations, preferred frames and the violation of the 
conservation of momentum. However, even in the Earth-Moon case, or the Solar 
system in general, the self-energies involved are still small and do not test the 
SEP in strong-field regimes where deviations to higher order terms of e could be 
present. It is at that point where the circular pulsar- white dwarf systems become 
important. 

For neutron stars, e ~ 0.15, which is large, in particular considering e — 0.5 
for a black hole, and much larger than the self-energy of a white dwarf, e ~ 10~ 4 . 
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Therefore, the pulsar and white dwarf companion of a binary system should feel 
a different acceleration due to the external Galactic gravitational field if SEP 
is violated. Similar to the Nordvedt effect, this should lead to a polarisation of 
the pulsar- white dwarf orbit ("gravitational Stark effect" 28 ). The eccentricity 
vector of such a binary system should therefore have two components, one con- 
stant component due to the external acceleration, and another one that evolves 
in time following relativistic periastron advance. Since the present direction of 
the evolving eccentricity vector is unknown, a careful analysis of all relevant sys- 
tems in a statistical manner is needed. Significant contributions to the results 
are made by long orbital-period and small-eccentricity systems, i.e. where Pb/e 2 
is large |77| . 

Preferred- Frames & Conservation Laws Some metric theories of gravity 
violate SEP specifically by predicting preferred-frame and preferred-location ef- 
fects. A preferred universal rest frame, presumably equivalent with that of the 
Cosmic Microwave Background (CMB), may exist if gravity is mediated (in 
part) by a long-range vector field. This violation of local Lorentz-invariance is 
described by the two PPN parameters a± and ct2- While both parameters can 
be tightly constrained in the weak-field limit of the solar system, ot\ can also be 
studied in the strong-field regime by analysing the same low-eccentricity pulsar- 

1/3 

white dwarf systems with a figure-of-merit given by P b /e: If a\ were different 
from zero, a binary system moving with respect to a preferred universal rest 
frame would again suffer a long-term change in its orbital eccentricity. In a sta- 
tistical analysis similar to that for the study of the gravitational Stark-effect, 
the binary 5.3-ms pulsar PSR J1012+5307 is particularly valuable. It not only 
has an extremely small orbital eccentricity, for which only an upper limit of 
e < 8 x 10 -7 (68% C.L.) was found from Jodrell Bank and Effelsberg observa- 
tions [29] . but its optically detected white dwarf companion also provides full 
3-d velocity information relative to the CMB. Using this and other systems, Wex 
P?l3"| derives \a x \ < 1.2 x 1(T 4 (95% C.L.) which is slightly better than the solar 
system limit (see [2]). 

In cases where theories both violate the Lorentz-invariance and the conser- 
vation of momentum, the equation of motion for a rotating body in the post- 
Newtonian limit contains so-called self-acceleration terms. This self-acceleration 
of the body's centre depends on the internal structure of the rotating body and 
results from the breakdown in conservation of total momentum. Another term 
in the self-acceleration involves the body's motion relative to a universal rest 
frame. Both contributions relate to the PPN parameter a 3 that can be tested 
using pulsars as isolated rotating objects |3UI31| . or as bodies in binary systems 
where both pulsar and companion suffer self-acceleration, leading to polarized 
orbits |32| . The limits derived by this second method using circular pulsar- white 
dwarf systems are much tighter than studying the spin periods of isolated pul- 
sars, resulting in a 3 | < 1.5 x 10~ 19 (95% C.L.) c 3f.. This result for 03 and a 
limit set on (0:3 + £2) < 4 x 10 -5 jSS] constrains the PPN parameter £2, which 
describes the non-conservation of momentum. Derived from a limit on the sec- 
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ond period derivative, P, of PSR B1913+16 and, hence, its acceleration, the 
interpretation of this limit may be complicated as non-zero Ps could arise from 
a number of sources such as timing noise fjj. 



Gravitational Dipole Radiation Essentially any metric theory of gravity 
that embodies Lorentz-invariance in its field equations predicts gravitational ra- 
diation. However, the details of these predictions may differ in the speed of the 
gravity waves, the polarization of the waves and/or the multi-polarities of the 
radiation. If a theory satisfies SEP, like general relativity, gravitational dipole 
radiation is not expected, and the quadrupole emission should be the lowest mul- 
tipole term. This arises because the dipole moment (centre of mass) of isolated 
systems is uniform in time due to the conservation of momentum and because 
the incrtial mass that determines the dipole moment is the same as the mass that 
generates gravitational waves. In alternative theories, while the inertial dipole 
moment may remain uniform, the gravity wave dipole moment may not, since in 
a violation of SEP the mass generating gravitational waves depends differently 
on the internal gravitational binding energy of each body than does the inertial 
mass pp. If dipole radiation is predicted, the magnitude of this effect depends 
on the difference in gravitational binding energies, expressed by the difference in 
coupling constants to a scalar gravitational field, (a p — a c ). For a white dwarf 
companion |a c | <C \a p \, so that the strongest emission should occur for short- 
orbital period pulsar-white dwarf systems. Again, the binary pulsar J1012+5307 
becomes extremely useful. Given its vanishing eccentricity, the change in orbital 
period due to dipole radiation becomes 

pdipo i e „ ^ 2 G* M p M c 2 AM 

^ " c*P b M p + M c ap +U \ ( *J , [ ] 

where G* is the "bare" gravitational constant. With the optically detected com- 
panion, the measured radial velocity can be use to correct for Doppler effects. 
For this system gjg, Wex derives a limit of |a p | 2 < 4 x 1(T 4 (95% C.L.). 



Time Variability of the Gravitational Constant Three different pulsar 
tests are available to test the time variability of the gravitational constant, G, 
and to derive an upper limit for G/G. A time variability is only allowed if SEP 
is violated due to preferred locations in space and time. In the case of pulsars, 
a changing G would change the gravitational binding energy of a neutron star 
and thereby possibly also its moment of inertia, which would cause a change 
in the spin-down behaviour, namely a contribution to P. A comparison with 
observed values leads to limits of the order of G/G < 10 -11 yr _1 4 . A slightly 
more stringent limit can be derived from the effects that a varying G would 
have on orbital periods |34I35| . Both limits are still about an order of magnitude 
above the limits set by solar system tests. Moreover, they depend to some degree 
on the compactness of the neutron star and its equation of state, so that they 
are not truly theory independent £Q. An interesting alternative test uses the 
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mass determination for neutron stars 36 , utilising that the Chandrasekhar mass 
depends directly on G. Studying the mass of millisecond pulsars as function of 
pulsar age, a strong limit of G/G < (-0.6±4.2)1CT 12 yr" 1 (95% C.L.) is derived. 
However, while the mass of neutron stars can be determined quite accurately in 
relativistic binaries (see Section l5~2*)l . an age estimation relying on Eq. J2J) can 
contain considerable uncertainty. 

5.2 Tests using Double Neutron Stars 

Even though in all metric theories matter and non-gravitational fields respond 
only to the space-time metric, it is possible that scalar or vector fields exist 
in addition to the metric. Damour & Esposito-Farese developed a framework to 
study theories at a second post-Newtonian (2PN) level where gravity is mediated 
by a tensor field and one or more scalar field |38| . These theories are interesting 
since scalar partners to gravitons arise naturally in quantum gravity and unified 
theories. Damour & Esposito-Farese show that it is possible to construct corre- 
sponding theories where deviations from general relativity are not visible in the 
weak field but only manifest themselves in a "spontaneous scalarization" if the 
strong field limit is approached. They conclude that current solar system tests 
and also upcoming satellite missions will not be able to replace the strong field 
tests provided by radio pulsars. Indeed, they use the DNSs, PSR B1534+12 and 
PSR B1913+16, together with PSR B0655+64 and solar system tests, to signif- 
icantly constrain the parameters describing the coupling of matter to the scalar 
field More stringent limits have been presented recently using also results 
for the pulsar-white dwarf system PSR J1141— 4565 40 . A more classical ap- 
proach using DNSs for tests of theories of gravity is made with the measurement 
of post-Keplerian (PK) parameters as observables. 

Because of the strong gravitational fields, we expect DNSs to suffer large 
relativistic effects. In this case, we cannot necessarily assume that we under- 
stand the underlying physics, even though general relativity appears to describe 
the physics in the solar system to high precision. Therefore, a purely theory- 
independent approach like the PPN approximation is difficult to realize. Instead, 
one can only use an existing theory of gravity and check if the observations are 
consistently described by the measured Keplerian and PK parameters. In each 
theory, for point masses with negligible spin contributions, the PK parameters 
should only be functions of the a priori unknown pulsar and companion mass, 
M p and M c , and the easily measurable Keplerian parameters. With the two 
masses as the only free parameters, an observation of two PK parameters will 
already determine the masses uniquely in the framework of the given theory. 
The measurement of a third or more PK parameters then provides a consistency 
check. In general relativity, the five most important PK parameters are given to 
lowest Post-Newtonian order by (e.g. 0): 




(15) 
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^2/3 (P b V /3 McjMp + 2M e ) 

• _ 1927T 5/3 /nr 5/3 (1 + I g2 + I e4 ) M P M C 
b 5 y2?r y (l-e 2 )V2 (Af„ + M c )V3' 

r =T M C , (18) 

where Ph is the period and e the eccentricity of the binary orbit. The masses 
Alp and M c of pulsar and companion, respectively, are expressed in solar masses 
(M ). We define the constant T Q = GM Q /c 3 = 4.925490947/^s where G de- 
notes the Newtonian constant of gravity and c the speed of light. The first PK 
parameter, cj, is the easiest to measure and describes the relativistic advance 
of periastron. According to Eq. (|15f) it provides an immediate measurement of 
the total mass of the system, (M p + M c ). The parameter 7 denotes the ampli- 
tude of delays in arrival times caused by the varying effects of the gravitational 
redshift and time dilation (second order Doppler) as the pulsar moves in its el- 
liptical orbit at varying distances from the companion and with varying speeds. 
The decay of the orbit due to gravitational wave damping is expressed by the 
change in orbital period, Pb- The other two parameters, r and s, are related 
to the Shapiro delay caused by the gravitational field of the companion. These 
parameters are only measurable, depending on timing precision, if the orbit is 
seen nearly edge-on. 

Until very recently, only two DNSs had more than two PK parameters de- 
termined, the 59-ms pulsar B1913+16 and the 38-ms PSR B1534+12. For PSR 
B1913+16 with an eccentric (e = 0.61) 7.8-hr orbit, the PK parameters Cj, 7 and 
Pb are measured very precisely. Correcting the observed Pb value for effects of 
relative motion (see Section 14.10 . the measured value is in excellent agreement 
with the prediction of general relativity for quadrupole emission (see Fig.0J). This 
result demonstrates impressively that general relativity provides a self-consistent 
and accurate description of the system which can be described as orbiting point 
masses, i.e. the structure of the neutron stars does not influence their orbital 
motion as expected from SEP. The precision of this test is limited by our knowl- 
edge of the Galactic gravitational potential and the corresponding correction 
to Pb- The timing results for PSR B1913+16 provide us with the most precise 
measurements of neutron star masses so far, i.e. M p = (1.4408 ± 0.0003)Mq 
and M c = (1.3873 ± 0.0003)M Q j^Z]. It is worth pointing out that these values 
include the unknown Doppler factor. 

The 10-hr orbit of the second DNS PSR B1534+14 (e = 0.27) is observed 
under fortunate circumstances, it is seen nearly edge-on. Thereby, in addition to 
the three PK parameters observed for PSR B1913+16, the Shapiro-delay param- 
eters r and s can be measured, enabling non-radiative aspects of gravitational 
theories to be tested, as Pb is not necessarily needed. In fact, the observed value 
of Pb seems to be heavily influenced by Shklovskii-terms, so that the correspond- 
ing line fails to meet the others in a M p -M c diagram. However, assuming that 
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Fig. 4. (left) Shift in the periastron passage of the DNS PSR B1913+16 plotted as 
a function of time, resulting from orbital energy loss due to the emission of gravita- 
tional radiation. The agreement between the data, now spanning almost 30 yr, and 
the predicted curve due to gravitational quadrupole wave emission is now better than 
0.5%. Figure provided by Joel Weisberg and Joe Taylor, (right) "Mass-mass" diagram 
showing the observational constraints on the masses of the neutron stars in the double- 
pulsar system J0737-3039. The shaded regions are those which are excluded by the 
Keplerian mass functions of the two pulsars. Further constraints are shown as pairs of 
lines enclosing permitted regions as predicted by general relativity: (a) the measure- 
ment of di gives the total system mass mA + niB = 2.59 Mq; (b) the measurement of the 
mass ratio R = ruiA/ms = 1.07; (c) the measurement of the gravitational redshift/time 
dilation parameter 7; (d) the measurement of the two Shapiro delay parameters r and 
s. Inset is an enlarged view of the small square encompassing the intersection of the 
three tightest constraints, representing the area allowed by general relativity and the 
present measurements. 

general relativity is the correct theory of gravitation, the deviation from the 
predicted value and the measured proper motion, fi, can be used to compute the 
necessary correction and hence the distance to the pulsar, d = 1.02 ± 0.05 kpc 

ED 

5.3 Tests Using Profile Structure Data 

In addition to the use of pulsars as clocks, strong gravity effects can also be 
tested using pulse structure data, namely the effects of "geodetic precession" in 
the DNSs PSR B1913+16 and PSR B1534+14. In both cases, the pulsar spin 
axis appears to be misaligned with the orbital angular momentum vector. In such 
a case, general relativity predicts a relativistic spin-orbit coupling, analogous 
to spin-orbit coupling in atomic physics. The pulsar spin precesses about the 
total angular momentum, changing the relative orientation of the pulsar towards 
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Earth. As a result, the angle between the pulsar spin axis and our line-of-sight 
changes with time, so that different portions of the emission beam are observed 
|42| . Consequently, changes in the measured pulse profile and its polarization 
are expected. In extreme cases, the precession may even move the beam out of 
our line-of-sight and the pulsar may disappear as predicted for PSR B1913+16 
for the year 2025 03]. See the review by Kramer f° r a detailed description 
of this effect and its observations. 

5.4 Recent Discoveries 

The tremendous success of recent surveys, in particular those using the Parkes 
telescope (e.g. @H]), has lead not only to the discovery of more than 700 new 
pulsars, but also to some very exciting new binary systems. Until recently, only 
five DNSs were known. This situation has changed and now eight systems can be 
studied. The most recent addition is a DNS discovered in the Parkes Multibeam 
(PM) Survey, PSR J1756-2251 (Faulkner et al, in prep.). This system shows 
similarities with PSR B1913+16 as its orbital period is somewhat less than 8 
hours although its eccentricity is smaller (e = 0.18). 

The recent discoveries benefit from larger available computing power which 
enables so-called "acceleration searches" for fast orbiting binary pulsars. Such 
techniques try to correct for the usually made assumption that the pulse pe- 
riod remains constant during the observations. For compact binary systems, this 
assumption is violated due to large Doppler shifts in period, resulting in much 
reduced sensitivity in standard Fourier searches. The employment of acceler- 
ation codes has therefore lead to a number of new binary pulsars with short 
orbital periods. Another example is PSR J1744— 3922 (Faulkner et el., in prep.). 
This 172-ms pulsar is in an almost circular 4.6-hr orbit and hence only the 
second long-period pulsar in such a short orbit. The other such pulsar is PSR 
Jl 141— 6545 which is a 393-ms PM pulsar in an eccentric 4.5-hr orbit Both 
pulsars appear to have a white dwarf companion, but while PSR J1141— 6545's 
companion is heavy [M c ~ IM@), the new pulsar's companion is probably much 
lighter (M c > 0.08M©). Whilst these are indeed exciting discoveries, the most 
stunning success is clearly the recent discovery of the first double-pulsar system, 
J0737-3039. 

6 The Double-Pulsar 

The 22.8-ms pulsar J0737-3039 was discovered in April 2003 021- It was soon 
found to be a member of the most extreme relativistic binary system ever dis- 
covered: its short orbital period (P& = 2.4 hrs) is combined with a remarkable 
high value of periastron advance (u> = 16.88 ± 0.09deg/yr, i.e. four times larger 
than for PSR B1913+16!) and a short coalescing time (<~ 85 Myr). The latter 
time-scale boosts the hopes for detecting a merger of two neutron stars with 
first-generation ground-based gravitational wave detectors by about an order of 
magnitude compared to previous estimates based on only the DNSs B1534+12 
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and B1913+16 07J- Consequently, during the lecture I had presented this pulsar 
already as the most beautiful laboratory for testing general relativity found so 
far, pointing also out that with a geodetic precession period of only 70 yr future 
studies should reveal interesting and exciting results. But little did we know then 
which surprise was still waiting for us. 

In October 2003, our team detected radio pulses from the second neutron star 
when data sets covering the full orbital period were analysed |48j . The reason why 
signals from the 2.8-s pulsar companion (now called PSR J0737— 3039B, hereafter 
"B") to the millisecond pulsar (now called PSR J0737-3039A, hereafter "A") 
had not been found earlier, became clear when it was realized that B was only 
visible clearly for two short parts of the orbits. For the remainder of the orbit, 
the pulsar B is extremely weak and only detectable with the most sensitive 
equipment. The detection of a young pulsar-companion B clearly confirmed the 
evolution scenario presented in Section [21 and made this already exciting system 
sensational, providing a truly unique testbed for relativistic gravity. 

Indeed, we have now measured A's Co and 7 and we have also detected the 
Shapiro delay in the pulse arrival times of A due to the gravitational field of B, 
providing a precise measurement of the orbital inclination of sini = 0.9995(^g 2 ). 
Obviously, as another strike of luck, we are observing the system almost com- 
pletely edge-on which allows us to also probe pulsar magnetospheres for the very 
first time by a background beacon. The measurements already provide four mea- 
sured PK parameters, resulting in a wia — tub plot shown in Fig.0] The orbital 
decay due to gravitational wave emission is already visible in the data, but the 
uncertainties are yet too large to provide a useful constraint. However, in addi- 
tion to tests with these PK parameters, the detection of B as a pulsar opens up 
opportunities that go well beyond what has been possible so far. With a measure- 
ment of the projected semi-major axes of the orbits of both A and B, we obtain a 
precise measurement of the mass ratio, R(ttia, tob) = mx/m^ — ieb/^a, provid- 
ing a further constraint displayed in Fig.0] For every realistic theory of gravity, 
we can expect the mass ratio, R, to follow this simple relation 02]. Most im- 
portantly, the i?-line is not only theory- independent, but also independent of 
strong-field (self-field) effects which is not the case for PK-parameters. This pro- 
vides a stringent and new constraint for tests of gravitational theories as any 
intersection of the PK-parameters must be located on the i?-line. At the same 
time, it provides us already with very accurate mass measurements for the neu- 
tron stars, M A = (1.337 ±O.OO5)M and M B = (1.250 ±O.OO5)M , respectively, 
making B the least-massive neutron star ever observed. 

The equations for the PK parameters given in Section 15.21 are all given to 
lowest Post-Newtonian order. However, higher-order corrections may become 
important if relativistic effects are large and timing precision is sufficiently high. 
Whilst this has not been the case in the past, the double pulsar system may 
allow measurements of these effects in the future 0H]- One such effect involves 
the prediction by general relativity that, in contrast to Newtonian physics, the 
neutron stars' spins affect their orbital motion via spin-orbit coupling. This effect 
would be visible clearest as a contribution to the observed wina secular 1501 and 
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periodic fashion For the J0737— 3039 system, the expected contribution is 
about an order of magnitude larger than for PSR B1913+16, i.e. 2 x 10 -4 deg 
yr _1 (for A, assuming a geometry as determined for PSR B1913+16 02]). As the 
exact value depends on the pulsars' moment of inertia, a potential measurement 
of this effect allows the moment of inertia of a neutron star to be determined for 
the first time |52| . Obviously, the double pulsar system offers improved but also 
new tests of general relativity. The current data already indicate an agreement 
of the observed with the expected Shapiro parameter of s bs/sexp = 1.00007 ± 
0.00220 (Kramer et al. in prep.) where the uncertainties are likely to decrease. 

7 Conclusions & Outlook 

Millisecond pulsars find a wide range of applications, in particular for precise 
tests of theories of gravity. After the discovery of pulsars thereby marked the 
beginning of a new era in fundamental physics, pulsars discovered and observed 
with the future Square-Kilometer-Array (SKA) will further transform our un- 
derstanding of gravitational physics. The SKA's sensitivity will discover the 
majority of pulsars in the Galaxy, almost certainly providing the discovery of 
the first pulsar-black hole system. For tests of general relativity such a system 
would have with a discriminating power that surpasses all its present and fore- 
seeable competitors [HHl- In particular, we could directly test BH properties as 
predicted by general relativity, such as the Cosmic Censorship Conjecture or the 
"no-hair" theorem. Moreover, the pulsars discovered with the SKA would act as 
arms of a huge gravitational wave detector enabling the study of a possible grav- 
itational wave background in a frequency range that is inaccessible to LIGO or 
even LISA. Clearly, the SKA will provide yet another leap in our understanding 
and application of pulsars. 
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